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ABSTRACT 

This letter extends our previous fluctuation analysis of the near infrared back- 
ground at 1.6/im to the 1.1/im (FllOW) image of the Hubble Ultra Deep field. 
When all detectable sources are removed the ratio of fluctuation power in the 
two images is consistent with the ratio expected for faint, 2; < 8, sources, and 
is inconsistent with the expected ratio for galaxies with z > We also use 
numerically redshifted model galaxy spectral energy distributions for 50 and 10 
million year old galaxies to predict the expected fluctuation power at 3.6/um and 
4.5^m to compare with recent Spitzer observations. The predicted fluctuation 
power for galaxies at z = 0-12 matches the observed Spitzer fluctuation power 
while the predicted power for z > 13 galaxies is much higher than the observed 
values. As was found in the 1.6/xm (F160W) analysis the fluctuation power in 
the source subtracted FllOW image is two orders of magnitude below the power 
in the image with all sources present. This leads to the conclusion that the 0.8- 
1.8|Um near infrared background is due to resolved galaxies in the redshift range 
z < 8, with the majority of power in the redshift range of 0.5-1.5. 

Subject headings: cosmology: observations - diffuse radiation - early universe, 
galaxies: high-redshift 
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Introduction 



In a previous paper (IThompson et al.l 120071 ) we addressed tlie nature of tlie Near In- 
fraRed Background (NIRB) at 1.6;um using NICMOS observations in the Hubble Ultra Deep 
Field (HUDF). The area covered by the NICMOS observations comprises the NICMOS Ultra 
Deep Field (NUDF) which is smaller than the HUDF. The NIRB at l.l^m and l.Gfim has 
a total power of 7 nWm~^sr~^ emitted by resolved galaxies, pred ominantly in the r e dshift 
range z = 0.5-1.5. This is in contrast to the previous results from iMatsumoto et al.l (120051 ) 
that claimed a peak flux of 70 nWm~^sr"^ at 1.4yum. The discrepancy arose not in the value 
of the total measured flux but in the attribution of the f lux to the components o f zodia- 
cal emission, resolved stars and galaxies and an excess. In IMatsumoto et al.l (120051 ) models 
were used to determine the first two components and the difference bet ween the modeled 



compo nents and the total observed flux was attributed to an excess. The IMatsumoto et al. 



( 120051 ) measurements did not extend shortward of 1.4^m and the sharp drop in the flux 
levels betwe en the 1.4^tm value and sho rter wavelengths was initially attributed by some 



authors (eg. ISalvaterra fc Ferraral (120031 )) as due to the Lyman break in very high redshift 



galaxies. Our analysis measured the flux from the resolved stars and galaxies and the zodia- 
cal background and found no excess. The difference was in the amount of flux attributed to 
the zodiacal background . Our measured value was higher than the modeled value used by 
Matsumoto et al.l (120051 ) and the difference was equal to the flux attributed to an excess. We 
therefore concluded that unless there was a extragalactic flux component that was extremely 
fiat, mimicking zodiacal flux, there is no excess and the near infrared background is resolved. 



Thompson et al.l (120071 ) also addressed the spatial fluctuations in the NUDF field relative 



to the fluctuatio n s fou n d in deep calibration 2M ASS images at the same wavelength by 
Kashlinsky et al.l (120021 ) . iKashlinsky et al.l ((2002) subtracted all of the resolved sources from 
the 2MASS images and found fluctuations in excess of that expected by shot noise at long 
wavelengths. The excess was attributed to fluctuation power from galaxies at red shifts 



of 10 and above in several studies (eg. iMagliocchetti. Salvaterra fc Ferraral (120031 )). A 



fluctuation analysis of the NUDF F160W image showed that when sources were removed 
from the NICMOS image down to the level that were removed froni the 2M ASS image, the 
fluctuations were consistent with those found by IKashlinsky et al.l (120021 ). When sources 
were removed down to the detection limit of the combined ACS and NICMOS observations 
in the NUDF the fluctuation power dropped by almost two orders of magnitude. Since all of 
the detected sour c es had redshifts between and 8, this meant that the fluctuations found by 
Kashlinskv et al.l (120021) wer e frorn galaxies in the redshift range between and 8. Further 
analysis in iThompson et al.l (120071 ) of images that only retained sources in a given redshift 
bin determined that the majority of fluctuation power came from galaxies in redshift range 
between 0.5 and 1.5. The fluctuation spectra from galaxies with redshifts greater than 4 
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were not detectable above the noise. 



More r ecent fluctuation analyses of Spitzer IRAC images at 3.6/im and 4.5/im (IKashlinsky et al. 



(I2OO5I ) and IKashlinsky et al.l (l2007bl )) interpret the residual fluctuations after source sub- 



tractio n as being due to high redshift, possibly Population III, galaxies (IKashlinsky et al. 



2007al ). This assertion appears to be based on the assumption that all galaxies below the 
IRAC detecti on limit must be fain t due to their high redshift. This interpretation has been 
challenged by ICooray et al.l (120071 ) who did further source subtraction in the fleld using deep 
ACS images to identify faint sources missed by the IRAC images. Th e y foun d the removal 
of these sources greatly reduced the signal seen by IKashlinsky et al.l (l2007bl ) and the sig- 
nal must be due to sour ces with redshifts less than 8 that are visible in the ACS images. 
Kashlinsky et al.l (l2007al ) replied that the exces sive amount o f area removed from the IRAC 
images in th e source subtraction car ried out by lCooray et al.l (|2007l ) invalidates their result. 
Modehng by ISalvaterra et al.l (120061 ) indicates that the observed IRAC fluctuations may be 
consistent with those created by Pop II stars at redshift greater than or equal to 5. The 
present study is motivated in part by a desire to investigate whether the observations in 
the NUDF at both 1.1/im and 1.6/xm are consistent with the high redshift interpretation of 
the the IRAC observations by examining the color of the fluctuations and their expected 
extrapolation to the IRAC bands. 



2. Data and Analysis 



The location, size and pixel scal e of th e FllOW image is identical to the F160W NUDF 
image analyzed i n iThompson et al.l (120071 ). Details of the image preparation are given in 
Thompson et al.l (120051 ). The basic image production after the processing of the individ- 
ual images is production of a background image which is the median of all of the individ- 
ual images, subtraction of the backgro und from the individual i mages and then combining 
the images with the drizzle procedure (jFruchter and Hookll2002l ). The fluctuation analysis 
of the FllOW image is identical to the analysis techniques used on the F160W image in 
Thompson et al.l (120071 ) and is not repeated in detail h ere. Source subtraction is accom- 
plished using the SExtractor (IBertin and ArnoutI Il996l ) (SE) pixel map where each pixel 
that is part of a source is given a value equal to the sour ce ID number an d all o ther pixels 
have a value of zero. The photometric redshifts derived in lThompson et al.l (120061 ) were used 
to identify the redshift of each source. Due to the much higher resolution of the NICMOS 
images relative to the IRAC images, even th e all source subtr a cted im age re tained 93% of 



i ts pix els. This means that the objections of IKashlinsky et al.l (|2007al ) to the ICooray et al. 



(j2007 ) analysis do not apply here. All source pixels were set to zero in the subtracted image. 
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There was no attempt to replace them with random noise. The small area of the subtracted 
sources means that this had no affect on the final fluctuation spectrum. Figure [1] shows 
the fluctuation power at 1.1/im and 1.6;um. To determine the redshift distribution of the 
l.l^m fluctuations the fluctuation analysis was also performed on the l.l/xm image with 
all sources except those in a given redshift bin removed. The 1.1/im fluctuation power ver- 
su s redshift distribu t ion is essentially identical to the 1.6/im distribution shown in Figure 6 
of [Thompson et al.l (120071 ). There is a peak near z=l and no discernible power above the 
background for z > 4. 



3. Fluctuation Colors 



To investigate the nature of the sources of the fluctuation power observed in the all 
source subtracted NICMOS and IRAC images we calculate the expected fluctuation power 
color as a function of redshift using galaxy Spectral Energy Distributions (SEDs). At any 
given angular scale the fluctuation power for a given wavelength band is directly proportional 
to the power of the galaxies (lylu) in that band. The power ratios of the bands at any 
redshift are easil y calcu lated from the numerically redshifted model SEDs. The analysis of 
Thompson et al.l (120061 ) used 7 primary template SEDs labeled 1-7 going from early to late 
galaxies. The first 5 are o bserved SEDs and SEDs 6 and 7 are calculated from the models of 
Bruzual fc CharlotI (120031 ). The properties of these models are given in Table [H Template 6 
is for a 50 million year old galaxy and Template 7 is fo r a 10 million year old galaxy with the 
lowest metallicity available in the iBruzual fc CharlotI (120031) models . Tern plate 7 is similar 
to a Pop. Ill SED but redder than a true Pop. Ill SED (see ISchaererl (120021 ) for the expected 
differences). Although a Pop. Ill SED may be bluer than template 7, the sharp increase in 
the 3.6/im to 1.6/xm fluctuation power at high redshift is due to the Lyman break passing 
through the 1.6/im and is therefore relatively independent of the SED to the red of the break. 



Thompson et al.l (120061 ) interpolated between the SEDs in increments of 0.1 to produce 
a final set of 61 different model SEDs. Figure |2] shows the distribution of the template 
numbers for all sources in the NUDF with redshifts greater than 2.0. The predominance of 
sources concentrate on templates 6 and 7, therefore, we will take them as the most probable 
SEDs for sources fainter than our source detection limit. We use these SEDs and the effective 
response functions in the NICMOS and IRAC filters to calculate the expected ratios of the 
power at redshifts between and 1 5. Absorption i n the Inter-Galactic Medium is modeled 



according to the treatment given in iMadau et al. I (119961 ) 



The horizontal lines in Figure |3] show the o bserved ratios of the fluct uati on power in the 
all source subtracted images using; this analysis, iKashlinsky et al.l (l2007bl ) and lThompson et al 
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( 20071) ■ Since the excess fluctuation power attributed to high redshift galaxies by lKashlinsky et al 



( l2007bl ) is at angular scales greater than 30" the horizontal lines show the ratios observed at 



100", the largest NICMOS angular scale. 

The observed 1.1/im to 1.6/im fluctuation power ratio is consistent with the lower red- 
shift template 6 predictions, but is incompatible with the calculated ratios at 2; > 8 for 
either of the templates. This is strong evidence that the source subtracted fluctuations in 
the NICMOS images are due to galaxies at 2 < 8. The observed 3.6/xm to l.G^um ratio lies 
between the ratios predicted by templates 6 and 7. Given the good match to the NICMOS 
flux ratios for template 6 we would expect higher 3.6^m fluctuations than observed. The 
difference may be due to cosmic variance since the flelds are different. It may also be that 
the much broader Spitzer PSFs subtract more background flux than the HST PSFs. For 
2; > 12 the predicted ratios become incompatible with observations. The 4.5/xm to 3.6/im 
ratio is essentially independent of redshift at the redshifts considered and therefore does not 
contribute any information on the redshift epoch of the fluctuations. The apparent agree- 
ment between the predicted and observed ratios for those bands, however, is further evidence 
that the observed fluctuations are due to faint galaxies at redshifts less than 8. 



3.1. Compatibility of the Data Sets 



An important question when comparing the NICMOS and IRAC source subtracted 
remnant fluctuations is wh e ther so urces are subtracted to the same level in both data sets. 
Altho ugh iKashlinskv et al" ( 2007b ) does not explicitly state the achieved source subtraction 
level, iKashlinsky et al.l (j2007al ) states that any remaining sources in the IRAC image must 
be fainter than 10-20 nJy. For the SEDs used in t his study the IRAC limi t corresponds to a 
NICMOS limit of 4.5-9 nJy. From Figure 3b from [Thompson et al.l (120071 ) the completeness 
limit is approximately 20 nJy at l.Gfim, which is 2 to 4 times brighter than the IRAC 
subtraction limit. The sources for subtraction, however, were identifled using detections 
in a combination of the ACS and NICMOS HUDF images, with the A CS images going 



significantly deeper due to their much longer integration time. Figure 7 in [Thompson et al. 



(120061 ) shows that the detection limit in the ACS F775W band is approximately 30 in AB 
magnitudes which is 3.7 nJy. We can use the SEDs discussed in fj3]to see what IRAC fiuxes 
these correspond to. The results are shown Figure |H Template SED 6 gives a limit of 20 nJy 
and the bluer template 7 SEP gives a limit of 10 ri Jy which are compatible with the IRAC 
subtraction limits given by lKashlinsky et al.l (l2007al ). The F775W band only removes sources 
at that limit for redshifts less than 3 but the ACS F850LP long pass filter removes sources at 
that limit for redshifts up to 6. The two data sets are compatible for all redshifts less than 
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6. At redshifts between 6 and 8 the NICMOS subtraction limit is not as stringent as the 
IRAC subtraction hmit. To account for this difference the horizontal solid line in Figures [3^ 
and b would have to be raised by a factor of 2 to 4 but that would still not be compatible 
with redshifts greater than 12. The redshift distribution of galaxies with NICMOS fluxes 
less than 20 nJy would have to be significantly different than galaxies above this limit for 
the z = 6-8 sources to contribute any power to the observed NICMOS fluctuations. 

3.2. Signal or Noise? 

All of the discussion above is predicated on the assumption that the fluctuations ob- 
served after source subtraction are due to faint galaxies below the detection limit and not 
noise sources from the instrument or telescope. To check for the signature of Gaussian dis- 
tributed noise in the NICMOS images we populated an image of the same size as the UDF 
FllOW and F160W images with random Gaussian distributed noise with a variance equal 
to that found in the real images. An identical fluctuation analysis gives the results shown 
by the triangles in Figure [TJ The random noise is close to the observed fluctuation power at 
small angular scales, but is far below the observed power at angular scales of 10" and larger, 
the region considered in this analysis. We therefore conclude that the observed fluctuation 
power is not due to random noise components such as read noise, Poisson photon statistics, 
or dark current in the NICMOS images. 

We originally did not include the FllOW image in the fluctuation analysis due to the 
possibility of residual flat fielding errors. The median of the 144 individual images measures 
any flat field errors very accurately. The FllOW mosaic image created by replacing all of 
the individual actual images with the median image does appear to contain a component of 
noise due to small flat field errors. It is evident as the higher fluctuation power in Figure [1^ 
at spatial scales between 10" and 80". A similar signature is not present in the l.l/xm source 
subtracted fluctuation power. This indicates that the subtraction of the median image from 
each actual image before the mosaic accurately removes any flat field errors in the regions that 
have no detectable sources. These are the regions that contribute to the all source subtracted 
fluctuation power. We therefore exclude flat fielding errors as a source of fluctuation power. 

If the source subtracted fluctuations are due to faint, low redshift, galaxies then the 
l.l//m and 1.6/im source subtracted images should be spatially correlated. A simple check for 
correlation is to subtract the FllOW image from the F160W to test whether the fluctuations 
are reduced as would be expected if they are correlated. Since the exact scaling is not known 
we subtracted the FllOW image with scalings between and 2 in 0.1 increments. The log 
weighted sum of the fluctuation powers for angular scales greater than 30" versus the scale 
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factor is shown in Figure [5^. There is a broad minimum for a scale factor of 0.6, near the 
observed ratio of F160W to FllOW fluctuation of 0.75 at 100". Figure [5b shows the nested 
fluctuation spectra. This shows that the subtraction of the scaled FllOW image reduces and 
alters the fluctuation spectrum at larger angular scales. It also indicates that the two images 
are partially correlated after the subtraction of all known sources as would be expected if 
fainter undetected sources were present in the images. The lower boundary in Figure [5)d is 
the upper limit on fluctuation power contributed by high redshift, l.l^um dropout, galaxies 
plus any remaining noise in the images. This limit is still above the Gaussian noise spectrum 
shown in Figured^ and b. 

At the largest scale of 100" the residual fluctuation power in the F160W minus FllOW 
image is 0.2 nWm~^sr~^. If the residual is due to sources with redshifts less than 10 
then the expected power in the IRAC 3.6/im band is 0.02-0.2 nWm~^sr~^ depending on 
which SEP temp l ate is u sed. This is consistent with the 0.1 nWm^^sr^^ power observed by 



Kashlinsky et al.l (l2007bl ) at the same spatial scale. If the residual power is due to galaxies at 
redshifts greater than 8 then the predicted power is in the range of 0.04 to a very large num- 
ber depending on the SED and redshift. The predicted power starts to exceed the observed 
flux at a redshift of 11 for template 6 and at a redshift of 12 for template 7. At a redshift of 
13.5 the predicted power exceeds the observed 3.6/im power by a factor of 10 for both SEDs. 
The residual F160W minus FllOW power to 3.6//m power ratio is therefore consistent with 
emission from galaxies with redshifts less than 8 but cannot definitively rule out that the 
observed 3.6/im fluctuation power comes from galaxies in the redshift range between 8 and 
13. Redshifts beyond 13, however, are not consistent with the ratio if the difference between 
the lower limit in Figure [5]d and the Gaussian noise is due to real sources. It should be noted, 
however, that the predominant source of fluctuation power in the observed sources comes 
from galaxies at redshifts between 0.5 and 1.5. Having the residual F160W minus FllOW 
fluctuation power coming from galaxies with redshifts above 8 would require a remarkable 
suppression of the luminosity function of the 0.5-1.5 redshift galaxies at low luminosities. 

A final, but not as compelling, piece of evidence that the residual fluctuations are due 
to galaxies is that the l.l/xm to 1.6/im fluctuation power color is similar to that expected 
from z < 8 galaxies. We cannot completely rule out an unknown source of spatial power, 
which would be the same in both images due to their similar observation pattern and data 
analysis. However, in the absence of any known source of such spatial power we conclude 
that the NICMOS source subtracted fluctuations are due to faint, 2; < 8, galaxies below the 
detection limit. These galaxies become evident in the fluctuation analysis because it utilizes 
93% of the pixels while the detected source power is contributed by only 7% of the image 
pixels. 
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Without far more knowledge of the IRAC instrument and data analysis procedures than 
we currently possess we cannot make the same arguments for the faint galaxy origin of the 
IRAC fluctuations. We can say that the IRAC fluctuations are at the levels predicted for 
faint, z < 8, galaxies from the NICMOS observations. Since all four bands have ratios 
consistent with a faint, z < 8, galaxies origin we assume that the fluctuations in all four 
bands have extragalactic origins. We note that even if the IRAC fluctuations were due to 
noise, it would still support our conclusions that the IRAC fluctuations are not due to very 
high redshift galaxies. 



4. Conclusions 

Our basic conclusion is that the source subtracted fluctuations observed in the two 
NICMOS bands and the IRAC bands at 3.6/um and 4.5/um are consistent with faint, z < 8, 
galaxies that are below the individual detection limit in the respective images. We also 
conclude that the observed ratio of l.lfim to 1.6;um fluctuations is incompatible with a 
significant fraction of power from z > 8 galaxies and that the ratio of the IRAC 3.6//m to 
NICMOS 1.6/im fluctuation power is incompatible for galaxies at z > 13. In light of these 
observations we find no evidence that any of the observed fluctuations require flux from 
high redshift galaxies. Of course such galaxies may exist, but we conclude that they do not 
contribute measurable power to the current observations. 



The argument in iKashlinsky et al.l (j2007al ) for the IRAC fluctuations being due to high 



redshift galaxies appears to be based on the assumption that any sources fainter than their 
detection limit must be at high redshift because they are faint. Examination of Figure 3b in 



Thompson et al.l (120071 ) shows that the z = 0-7 galaxy luminosity function shows no break 
up to a completeness limit of 5 nJy in the F160W band. This implies that there are many 
more low redshift galaxies fainter than the detection limit, i.e. galaxies that are not faint 
because they are at z > 8. 

A caveat relevant to this study is that if it is postulated that the IRAC fluctuations 
are due to galaxies with z > 13 then the predicted 1.6/im fluctuations would be below our 
detection level due to the very high ratio of IRAC power to NICMOS power at high redshifts. 
However if these high-redshift galaxies are the only source of the IRAC fluctuations, then 
a mechanism must be found to suppress the low redshift IRAC fluctuations predicted from 
the NICMOS residual fluctuations. At this time we know of no such mechanism. 



This article is based on data from observations with the NASA/ESA Hubble Space 
Telescope obtained at the Space Telescope Science Institute, which is operated by the Asso- 
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Fig. 1. — The fluctuation spectra from a) tlie FllOW NUDF image (left) and b) the F160W 
image (right). The open squares are for no sources subtracted, the open diamonds are for all 
sources with a F160W AB magnitude of 20 or brighter subtracted. The asterisks are for all 
sources subtracted and the crosses are for the median background field described in the text. 
The open triangles are for a Gaussian distributed noise image with a standard deviation 
equal to that of the images. The da shed line indicates the le vel of 1.6/im fluctuations found 
in deep calibration 2MASS fields by iKashlinsky et al.l (120021 ). 
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Fig. 2. — The histogram of SED template numbers for all sources with redshifts greater than 
2.0. The interpolated templates go in increments of 0.1 between the primary template SEDs. 
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Fig. 3. — The predicted ratios of the powers, vly^ in the FllOW band to the F160W band, 
dashed hne; the the 3.6^m Spitzer band to F160W band, sohd hne; and the Spitzer A.bfim 
to 3.6/im band, dash dot hne; versus redshift. The straight horizontal hnes indicate the 
observed values of the ratios at an angular scale of 100" with the same line style code, a) 
(left) is for templateG and b) (right) is for template 7, the hottest SED, which is consistent 
with a possible Population III galaxy. 
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Fig. 4. — The figure shows the IRAC fluxes in the 3.6/im and A.bfim bands that correspond to 
a subtraction hmit of 3.7 nJy in the ACS F775W band as a function of wave length for the two 
templa tes (6 and 7) used in this study. The quoted subtraction level from iKashlinsky et al. 
(l2007bl ) is 10 to 20 nJy, consistent with subtraction levels in the figure. The sharp rise at 
approximately z=3 is when the Lyman break starts to enter the F775W band. 
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Fig. 5. — a) The left panels shows the ln(/c) weighted sum of fluctuation power for angular 
scales greater than 30" versus the scahng factor of the FllOW image subtracted from the 
F160W image. The minimum is at a scahng factor of 0.6. b) The right panel shows the 
nested fluctuation power spectra for the different scaling factors. The lower boundary is the 
limit for the contribution of high redshift sources plus other sources of spatial noise. 
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Table 1. Properties of the template SEDs. Both SEDs are taken from mod els of 



Bruzual fc CharlotI (j2003r ) based on the P adova 1994 models. The IMFs are from IChabrier 

J2003h . 



Template Number 


Metallicity 


IMF 


Low Mass Limit 


Upper Mass Limit 


6 


0.004 


Chabricr 


0.1 


100 


7 


0.0001 


Chabrior 


0.1 


100 



